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Abstract 

The low-energy behavior of the astrophysical S-factor for El direct radiative 
captures a(p,7)b leading to loosely bound final states (b=a+p) is investi- 
gated. We derive a first-order integral representation for S(E) and focus on 
the properties around zero energy. We show that it is the competition be- 
tween various effects, namely the remnant Coulomb barrier, the initial and 
final centrifugal barriers and the binding energy, that defines the behavior of 
the S{E 0). Contrary to previous findings, we prove that S{E — > 0) is 
not determined by the pole corresponding to the bound state. The derivative 
5'(0) increases with the increase of the centrifugal barrier, while it decreases 
with the charge of the target. For /, = // + 1 the increase of the binding 
energy of the final nucleus increases the derivative S"(0) while for li = If — 1 
the opposite effect is found. We make use of our findings to explain the low 
energy behavior of the S-factors related to some notorious capture reactions: 

^Be(p,7)SB, ^^N(p, 7)1^0, i^0(p,7)i^F, '^^Ne{p,-/f^Na and '^^Mg{p,'yf^Al. 
PACS Numbers: 25.40.Lw, 25.40.Ny, 27.20.+n 
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I. INTRODUCTION 



Direct radiative capture cross sections of charged particles a(p,7)b drop so sharply for 
the low energy region of stellar environment that it is often impossible to measure these 
rates at the appropriate energies in the laboratory Then, typically, the direct capture 
cross sections are measured at higher energies followed by an extrapolation down to zero 
energy. Therefore, the knowledge of the energy behavior of the cross section as ii^ — >^ is of 
crucial importance. 

When an incident charged particle approaches a target it should penetrate through the 
Coulomb and centrifugal barriers in the initial state. The capture to a loosely bound state 
is the transition from the initial scattering state to the tail of the bound state wave function 
in the final state. This tail is affected by the Coulomb and centrifugal barriers in the final 
state. Thus, the radiative capture process is affected by four barriers, two Coulomb and two 
centrifugal barriers, implicitly present through the four parameters ?7j, ?7j, /j, /j. All barriers 
have a similar effect on the radiative cross section: they decrease the probability of the 
reaction at sub-barrier energies. 

The cross section is usually factorized into the Gamow penetration factor, and the as- 
trophysical S-f actor: 

a{E) = '-^S{E). (1) 

The Gamow penetration factor is no more than the probability for an s-wave proton to 
penetrate through the pure Coulomb barrier (i.e. the Coulomb potential is extended down 
to the origin r = assuming a pointlike target). Naturally, this factor defines only the gross 
energy behavior. In fact, as the capture to loosely bound states occurs under the Coulomb 
barrier, at a large distance from the target, the effective penetration factor is larger than the 
Gamow penetration factor in Eq. (|l]) p[. We call remnant Coulomb barrier the remaining 
energy dependence of the astrophysical factor S(E) due to the initial Coulomb barrier, after 
removing the Gamow factor. This remnant effect behaves in an opposite way to the normal 
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barrier behavior, namely, it reduces the S-factor as the energy increases, to compensate the 
overestimation by the Gamow factor. 

It turns out that, depending on the system, S-factors can feature completely opposite 
behaviors as one approaches zero energy. Despite the long history of S-factor calculations 
and the numerous papers published on this subject, no satisfactory physical explanation for 
the different behaviors has been presented. In recent papers some aspects of the low 
energy behavior of the direct capture S-factors were investigated. In the low energy rise 
of the S-factor was attributed to the pole in the energy plane, located at E = —e (with e 
being the binding energy of the final nucleus b = a + p relative to proton threshold). In 
the integral form for the first three terms of the Taylor expansion were derived, using 
a potential model. In either work, the physical reasons for specific patterns in the S-factor 
energy behavior were not considered. 

If the rise of S{E — > 0) were in fact caused by the pole how could one understand the 
decreasing S{E 0) for ^^0(p, 7)^^F(gr. st.), a case where the pole is extremely close to 
the threshold? Why should the S-factor for ^^0(p, 7)^''F(2si/2) have a steeper slope near 
the threshold than the S-factor for ''Be(p, 7)^6? What are the crucial parameters defining 
the behavior of the S-factor for the direct captures to loosely bound states? Is it possible to 
define general thumb rules? It is clear that accurate numerical values for direct capture S- 
factors are trivially calculated with present day technology. For that reason analytical works 
have not been highly appreciated. Regardless, this work is a search for the understanding 
of the general physical features, based on analytical considerations. 

The vast majority of the direct proton radiative capture reactions with interest to As- 
trophysics, are dominated by -El transitions. We will thus concentrate on the El transition 
for the a + p 6 + 7 process. We also confine ourselves to loosely bound final states. 
The interest on the behavior of the S-factors for the direct captures reactions goes beyond 
p-shell nuclei. The heavier, sc?-shell nuclei ||^, when near the proton drip line and/or near 
shell closure, have very small Q-values for the radiative capture. Thus the S-factor is domi- 
nated by single resonant and direct captures. Consequently, the findings of this work have 
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implications for a great number of today's interesting astrophysical cases. 

In this paper we use an analytical/integral representation of the transition matrix element 
for the El process and impose the asymptotic approximation. Following, we identify the 
physical ingredients that affect the low-energy behavior of capture reactions to loosely bound 
states. In Section |I| the effect of the pole of S{E) at negative energy is specifically addressed. 
The threshold behavior is studied in Section |ITli The competition between the remnant 
initial Coulomb barrier, the final Coulomb, the initial and final centrifugal barriers, as well 
as the trace of the singularity at E = —e are detailed and discussed. The S-factor dependence 



on the target charge is also studied. Section ^ consists of an application to some important 
astrophysical capture processes. Finally, in section 0, a summary of the results is presented. 



II. THE S-FACTOR AROUND THE POLE 

For El captures to loosely bound states, the reactions under consideration become purely 
peripheral and a simple sohd sphere model is adequate. Such an approach is also used when 
calculating the non-resonant capture in the R- matrix method 0. The S-factor for the 
capture from an initial state of a + p with orbital angular momentum /j, to a final state Jf 
of b = a + p, with the angular orbital momentum If, can be written as 

S{E)=A'^{k,ror^'Pi^ 



X 



r, 



_ / drr^l^_,^,,,^+i/2(2/^r)[F,(fc,r)G',(fc,ro)-F,(A;,ro)G,(fc,r)]p, (2) 




where ^ is a factor depending on spin-angular characteristics, masses, charges and constants: 



/ I / I I \ / I / I I \ I r II 1 fyy^Lj 
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(2 J, + 1) (2 Jp + 1) 2 + 1 ' ' m, 



(L + l)(2L+2)^^_ , / 
^ L [(2L-M)!!P^ ^^'^"''^ """^ hc^- 

Here C is the asymptotic normalization coefficient of the tail of the projection of the bound 
state wave function of the final nucleus 6 onto the two-body channel a+p. The Sommerfeld 
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parameters for initial and final states, rji = ZaC^fiap/k and r]f = Zac'^ fiap/ k, are defined 
in terms of the initial relative momentum of a + p, = yj2 jiap E, and the bound state 
momentum, k = ^2 fiap£, respectively. Note that the initial relative energy is E and the 
energy of the final state is —e, so the photon has energy E + e, momentum = [E + e)/hc 
and multipolarity L. The penetration factor is defined in terms of the regular and irregular 
Coulomb functions at r = tq: 

Pi^ = kro/[E^Xk,ro)+Gl{k,ro)] (4) 

and W-rj f ,if +i/2i'^ i^r) is the Whittaker function describing the radial behavior of the bound 
state wave function at r > ro, where ro is the solid sphere radius. As usual, Za is the charge 
of a and flap is the reduced mass of the a + p system. The constants used are Xp = 0.2118 
fm for the proton Compton wave number and mc^ = 931.5 MeV for the mass atomic unit. 
The factor 10^ is introduced to provide the S-factor in keV b. 

When the binding energy of final state is small and the a—p relative energy is low, nuclear 
re-scattering a — p can be neglected. This corresponds to neglecting the term containing 
the irregular Coulomb function, Gf.{k,rQ). One should keep in mind that whenever precise 
values for S{E > 0) are needed, phase shifts should be taken into account. However, here 
we are only concerned with the qualitative behavior of S(E) around zero. In section V 
comments on the accuracy of our results are presented. Eq. then simplifies to: 

S{E) = A^ikro) {k^roY'^^'l^, Jdr r^ iy_,^,z,+i/2(2 r) F,(A;, r)p. (5) 

Both, Eq. (0) and Eq. (|]), can be used for the sake of our following arguments. 

Next, we realize that when A; — > 0, extending the integral in Eq. (|^) from ro to zero is 
a reasonable approximation, as then the dominant contribution comes from the asymptotic 
region. Replacing the bound state wave function by the Whittaker function down to zero 
radius is certainly not accurate for k > 0, but we have verified that the qualitative features 
for S{E) remain. 

Using this approximation and the asymptotic expansion of the Whittaker function, we 
can obtain an analytical result for S(E), which is not intended to be quantitative but will 
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contain the main energy dependence as ^ 0. In the following pages we perform the 
detailed derivation, but advise to follow onto Eq. (|T^) in case the reader is searching for a 
final result only. 

First, we use the asymptotic expansion of the Whittaker function (at r — > oo) [0: 

w ,„./,(2.r)-.^ {^(Il^3!Ml±M±j}l\ (6) 
with (/)o = 1 and = 1 ■ Z ■ (/ — 1) ■ ...(/ — (j — 1)). Substituting in Eq. @ we arrive at: 

I f^ {h-m)Alf + m+3)^ / drr^e-^ ^-^F,{k,r)\f . (7) 



Before the integration over r is performed, by including the explicit behavior of the regular 
Coulomb function, we can estimate the analytical behavior of the integral in each term of 
the asymptotic expansion, 

oo _ 

T = / drr^ —^—^—Fi.ik.r). (8) 
The regular Coulomb function is given by: 

FiXK r) = CiXVi) {krf+^ e''^'' ,F,[k + 1 - ir],, 2k + 2;2 tkr], (9) 

where i-Fi[/j + 1 — irji, 2/j + 2;2ikr] is a confluent hypergeometric function (see Appendix, 
Eq. (Q), and 

C,,,0.2..e-.W.E|±i^. (10) 

The asymptotic behavior of the Coulomb function as r — oo is well known: 

Fi^{k,r) ~ sin{kr - -k - r]i ln2kr + ai^), (11) 

with ai- = arg r(/j + 1 + zr/j). Although for the loosely bound states we need to know the 
Coulomb wave function only for r > tq, the physical meaning of the Ci.{Vi) coefficient comes 
from the behavior of Fi.{k, r) at small r: 
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Fjfc,r) ~Vjr/,)(A;r)^ (12) 

It is clear from Eq. (0) that Ci-{Vi) defines the probabihty of finding a charged particle in 
the vicinity of r = 0, under the influence of the Coulomb and centrifugal potentials. Note 
that 

OlM - 2-. . ^ ^pJ^ m + «^). (13) 

In brief, we have replaced the lower integration limit tq by zero, we have used the 
asymptotic expansion for the bound state Eq. @ and the explicit representation of the 
Coulomb scattering wave in Eq. (H). After performing the integration, the S-factor simplifies 
to: 

xlE^i^M-^ 2i^i[/. + l-^r7.,2-r/^ + L + /,-j,2/, + 2,— (14) 
where the expansion coefficients are given by: 

= ih - Vf)j ih + Vf+ J)j —r, ■ (15) 



We note that the expression for S{E) in Eq. (|1^) is completely general. In particular, 
the expression for the neutron overlap integral can be obtained by calculating the limit 
Vf, Vi 0. 



Substituting Eq. (|1^) in Eq. (|1^) one arrives at: 
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^ (2/t)2^/ |(fc-iK)2-^/+^+'»|2 '^t^Q ' \ 2k J 

2k 

X ^F,[li + l-ir]i,2-r]f + L + k-j, 2/, + 2; -f. (16) 

rh Z Kj 

This equation bears all the necessary features to investigate the behavior of the S-factor 
approaching E ^ —e. 



Beforehand, we make some obvious remarks concerning the pole in S(E). Taking into 
account the asymptotic behavior of Fi.{k,r) given in Eq. (|1T]), we conclude that the singu- 
larity of Tj on the physical k-half-plane {Imk > 0) is due to the divergence of the radial 
integral on the upper limit (r oo) at k ^ i k. More explicitly, Tj behaves as: 



T, oc / drr^— — e-^'^V'?' . (17) 



When k IK, then irjf —>■ i]f. Then, it is easy to verify that: 

Tj oc{k-t ny-^-^ for j <L + 1 (18) 
oc ln(A; — in) for j = L + 1 . (19) 

Note that Tj is finite at k = i k for j > L + 1. Besides the square of the integral, S(E) 
contains also a factor k'^^'^^ (see Eq.|^. Consequently it is the term j = in the expansion 
that generates the pole in S{E), as first indicated in 0. All higher order terms, j > 0, go to 
zero as k i K, . Let Aj{E) be the agglomerate of all energy dependent factors in Eq. (|^: 

A,{E) = ie^-^^ (fcro)(A;,ro)2^+i |T,f . (20) 

Then its behavior as k ^ in (or E —e) is: 

AAE) oc ,^ . for j <L + l (21) 

oc [ln(E + e)f [E + ef^^^ for j = L + 1 . (22) 

The term j=0 of Eq. ( pOD does not depend on /j. Hence, if this pole really did govern 
the behavior of the S-factor at zero energy, then SiE) should increase rapidly as — 0^, 
independently of We know from experiment that this is not the case. In fact, we will 
show that the rise of S(Ei) at — > 0^ is dictated by higher order terms while the first order 
term, the only one with a pole at E = — e, actually decreases as E ^ 0^. To this end we 
calculate the behavior of SiE) around threshold for positive energies, — > 0"*", but also for 
negative energies, all the way down to E ^ —e^. 

This equation bears all the necessary features to investigate the behavior of the S-factor 
approaching E —e. Let us begin with the term j=0: 
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1 2A; 
X I (^3— p^TTITI- + 1 - 2 - r^, + L + 2k + 2; ^3—]!^. (23) 

In order to find the analytical properties of the astrophysical factor around the pole, E —e, 
it is useful to perform an argument transformation relating the hypergeometric functions 
with arguments z and 1/z, before taking the relevant limit k —>■ in [0: 

2k h ■ 

2Fi[h + l-ir]i,2~r]f + L + k,2k + 2; -] 

k — t K 

ci{k-t Kf+^-i^i^ + C2 (A; - i Kf-'^f+^+^\ (24) 
Using this result in Eq. (p3D we conclude that the leading singular term foT k ^ i k, is: 



whichever k, If or L, consistent with the conclusions drawn from Eq. pO|). 

It is important to realize that the behavior of the hypergeometric function 2-^1 [k + 1 — 
irji, 2 — rif + L + li, 2/j + 2; z = 2k/ {k — i k)] is very different around the pole {z 00) 
compared to its behavior around threshold {z ^ 0). The analytical expression of S{E) at 
threshold will been derived in Section |T|. 

2^1 [/« + 1 — irji, 2 — rjf + L + li, 2li + 2; 2k/ {k — i k)] is an analytic function in the finite 
fc-plane (except for k = Hence, assuming that the cuts drawn from branching points 

k = ±i K, go to infinity, we can make an analytic continuation of 2-^1 [k + l — irji, 2 — r]f + L + 
li, 2li + 2; 2k/{k — i k,)] into the complex fc-plane. Using Eq. (p3|) we arrive at an expression 
for the S-factor in the vicinity of the pole. Note that the continuation of S(E) will not 
necessarily be analytical at A; = 0, as the modulus operation may introduce a discontinuity. 
For k = ip and E = —p'^/2 ^ap, Eq. ( pSf ) takes the form 



K 



[r(2/, + 2)P ' ' {2Kf^f 

^ YrS2Fi[k + l-r{„2-'nf + L + k,2li + 2,^^]f, (26) 



where t][ = Zac'^ Hap/p- 

In Fig.|l] we present the energy behavior of 5(0) (-E) for a{p,'y)b, for masses Aa = 7 and 
Ab = 8 and angular momentum /j = 0, // = 1 and L = 1. It mimics ^Be(p, 7)^B, but a 
different set of proton binding energies are used. The energy interval is —e < E < e. All 
5(0) (-E) are normalized to unity at zero energy. It is obvious that, even for the smallest 
binding energy, in which the pole is closest to threshold, S'(o)(-E') decreases as E ^ 0~^. 
Whenever there is a rise of S(E) around zero, it can only be due to the sum effect of higher 
order terms in Eq. (p!6D, and not the pole. 

III. THRESHOLD BEHAVIOR OF THE S-FACTOR 

In this section, we look at the dependence of the derivative of S{E) on /j and If, and 
analyze the competition between the remnant initial Coulomb barrier and the trace of the 
singularity E = —e, discussed in the previous section. We also investigate the dependence 
of S{E) on charge and binding energy. 

Using the integral representation of the initial scattering state, and the final bound state, 
in the asymptotic approximation, one can deduce an integral expression for the S-factor close 
to threshold. This derivation is presented in detail in Appendix ^ and the result we obtain 
is: 



Vf jJo 



X 







(27) 



e-o de'" [(k (1 + 2 s) + e)2 + A;2]'«+i 

where an energy independent constant B has been introduced (Eq. [A.14D , S{s) is the spectral 
function defined in Eq. ( [A.8D and m = If + L + 1 — li. 

We mentioned that, by taking only the first term of the expansion of the Whittaker 
function defined in Eq. (|]) we can arrive at a similar expression for ^(o) where S{s) should be 
replaced by So{s) defined in Eq. ( |A.15| ). Although 5(0) has a pole at k = i k, a.s k —* +0 



its first derivative S'(^q^{E) > for all binding energies, as can be verified in Fig.|l|. The only 
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difference between S{E) and S(^q){E) is in the spectral functions for W^rf^. ^1^+1/2(2 Kr), Eq. 



( |A.1| ), and ^^+1/2(2 Eq. ( [A. 21) . The spectral functions at s — > 00 are, from Eq. 
( CT and Eq. ( W^ , 

Sis) = [sil + s)yf [-^T'' ™ s^'f and 5(o)(s) = s'f+^f . (28) 
1 + s 

There are only two competing factors here: the centrifugal barrier (If) and the Coulomb 
barrier (rif) of the final-state. Three situations may occur: 

i) for // = rif, the S-factor is represented exactly by the first term: S{E) = S(o){E); 

ii) for If < rjf, 5(o)(s)/5(s) 00: the contribution of the large s values is more pronounced 
for 5(0) than for S{E); 

iii) for If > rjf, iS(s)/5(o)(s) ''-^^ 00: the contribution of the large s values is more pro- 
nouncened for S{E). 

In brief, if the final-state Coulomb barrier prevails, 5'(o)'(0) > 5"(0) otherwise, if the final- 
state centrifugal barrier prevails, 5"(0) > ^(^(^(O). 

Only two cases need to be considered for dipole transitions L = 1: m = lif/j = /j + l 
and m = 3 if /j = // — !. We consider them separately. 



A. m = 1 when li = If -\- 1. 

After taking the required derivative in Eq. (^), the expression for S{E) is 

U2L+1 k 

S{E)=ABr^f^- + '^') (29) 



^ ,=1 



with 



00 



« ' -J ' f2is,k) 

fi{s,E) 



Let us first make some general qualitative remarks concerning the behavior of the in- 
tegrand near threshold E 0^. These will be useful later. In all cases considered, the 
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integrand peaks at Sm > 0. This peak occurs for a larger Sm, for l + l I than for / / — 1, 
reflecting the effect of the centrifugal barriers in the initial and final state. In addition, the 
peak of the integrand also shifts to larger s-values when the energy increases. 

Remarkably, the energy dependent part of the S-factor depends on E rather than on k 
and can be written as a positive definite function: 



The part in square brackets decreases monotonically at E ^ (positive derivative) while 
the integral /(fc^) increases as — > (negative derivative). The sign of the first derivative 
of S{E 0^), or in other words, the low-energy behavior of the S-factor, is defined by the 
sign of the derivative of E{E). 

Next, we focus our attention on /(A;^). The integrand of /(fc^) is a positive definite 
function of s and k and the two functions in the integrand of Eq. (|30|), fi{k'^, s) and 
/2(s, k), decrease monotonically with increasing energy. Now the interplay of the various 
ingredients becomes more transparent. Clearly, /2(s, k"^) is a remnant of the initial Coulomb 
barrier, as well as the second term of function /i(s, k"^), rif/[{l + 2 s)^ + E/e]'"'^'^. 

Note that, as /j = // + 1, the initial orbital cannot be an s-wave. As increasing /», shifts 
the peak of the integrand to higher values of s, the integrand, and, hence, lik"^), becomes 
less sensitive to energy variations. If I'{E) < increases as /j = // + 1 increases, the factor 
^7^"*"^ itjLiW^ k"^ + '^^) eventually prevails leading to the positive derivative S\E). 

In addition, increasing e shifts the pole away from threshold, such that I{k'^) becomes 
less varying with energy around threshold. So, increasing the binding energy, decreases 



We can derive an expression for the S-factor near threshold factorizing the energy de- 
pendence. First we rewrite /(/c^) as 



F{E)= 



(31) 



J'(fc2 = 0)|. 



1 



(32) 



{E + e) 



I{Q) + I'{^)E + 0{[E/ef) 



(33) 
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where /(O) ^ 0. 

Let us estimate the impact of the pole singularity on /(/c^) in the physical region at 
E 0^. In the case of a pure second order pole, I{k'^) = const, —e /'(0)//(0) = 2 (the 
order of the pole). Otherwise, the Taylor expansion holds: /(A;^) = /(0)+/'(0) -E'+0((-E'/£:)^), 
where /'(0)//(0) > 0. Then we get -e /'(0)//(0) = [2 - e /'(0)//(0)]. If [(e /'(0)//(0)] > 0, 
which corresponds to our case, the power of the singularity in I{k'^) as 0^ is weaker 

than 2, i.e. £/'(0)//(0) > -2. 

The explicitly terms appearing in Eq. (|33D can be derived: 

OO 







/'(O) = — J ds [s{l + s)ff 





l + s' 3(l + 2s)2'>+6 
[3 (/, + 2) (1 + 2 s) + 2 r]f] e^^^f^^^+^'l (34) 



^ 

Finally, we can also write down the threshold behavior of S{E), as S{E) = S{0) + 
S'{0)E + O{[E/6f) where: 

and 



1 + 4 + .. + /2 P/'O 



(36) 



From Eq. (p6|), it is clear that the sign of the derivative of the S-factor at ^ is defined 



by the sign of 



3 + (l + 4 + .. + /2)/^2_2^|j,(0)//(0)| 



. Note that the larger is /j, the larger 



becomes the positive term in this factor. 

We now illustrate these ideas with a few examples which relate to physical cases. 

(i) h=l, lf = 0. 

In Fig.0 we present the behavior of S{E) for the /j = 1 — > // = capture ''Be(p, 7)"^B", 
spanning a few values for the binding energy of "^B". All S{E) are normalized to unity at 
E = 0. The smaller the binding energy, the larger the increase of S{E — > 0). The centrifugal 
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barrier in the initial state is not strong enough to overcome the negative derivative caused 
by the remnant Coulomb {rji) and the singularity a.t E = —e. In addition there is the 
dependence on the spectral function [j^]''-'', where rjf depends of the inverse square root 
of the binding energy. For e = 0.06, 0.137, 0.35, 0.6 MeV the Coulomb parameters are 
rjf = 2.4, 1.6, 1.0, 0.76, respectively. Replacing // = in Eq. p^ ) we see that, for larger rjf, 
a significant contribution to the integral comes from smaller s- values, reducing e /'(0)//(0) < 



and l/rif"^, so that 



3 + l/r]} + 2er{0)/I{0) 



becomes negative. On the other hand, for 



smaller rjf, which corresponds to the pole at E = —e moving further away from threshold, 
the slope £/'(0)//(0) and l/rjf'^ increase and, correspondingly, 3 + l/r/j + 2£:/'(0)//(0) 
can become positive. 

We also compare S{E) with 5(0) obtained retaining only the first term of the Whit- 
taker expansion (Fig.^). Contrary to S'{E) we find that all calculated cases have positive 
'^(o)(0^)- This confirms that the negative derivative of S{E 0) is not due to the nearby 
pole as concluded in Section y. As rjf > If = 0, 5'(o)'(0) is always larger than S"(0). 

(ii) h = 2,lf = 1. 

In Fig. ^ are the S-factors for the U = 2 Ij = 1 capture ^Be(p, 7)"®B" , assuming a set 
of binding energies for '"^B". Strikingly, all binding energies produce a positive derivative 
for 5(0). In this case the centrifugal barriers in the initial and final states overcome the 
remnant Coulomb barrier and the pole ai E = —e. Thus, even at the lowest binding 
energy, e = 0.06 MeV, 5"(0) is positive. When the binding energy increases, the pole 
moves away from threshold (decreasing e I'{0)/ 1{0)) and rjf decreases, which means that 
3 + S/r/j + 2£:/'(0)//(0) becomes larger, as well as 5"(0). The general behavior for Zj > 
2 ^ // > 1 is identical. 

Fig. ^ also shows the behavior of 5(0) (-E), represented by the grey lines. This figure is 
a clear illustration of the points mentioned in the beginning of this section. For the cases 
e = 0.06 and e = 0.137 MeV, r// > //(= 1) and then 5(o)'(0) > 5'(0). The binding e = 0.35 
MeV corresponds to r/j = 1 for which 5(o)(0) = 5(0). Finally, for e = 0.6 MeV, rjj < If and 
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consequently 5(o)'(0) < ^'(0). 



B. m = 3 when 1^ = If — 1. 



Directly from Eq. p7|) we get 



X 



rl3 1 „ , I. 2 

(37) 



lim I ds Sis) r—. A-.^ r..T^ 



^-oy ^ ^ de3 [(k[1 + 2s] +e)2 + A;2]'»+i 

The explicit expression for the third derivative of the integrand is tedious, but the analysis 
is in all similar to the case with m = 1 . The energy dependence of the S-factor has the same 
form as Eq. (|3lD , 

h 



F{E) 



j=0 



m'), (38) 



only the integral changes due to the third derivative, which we here represent as I^{k'^). We 
keep in mind that the initial state now can be an s-wave {li = 0). 5(0) and 5"(0) are given 
by expressions identical to Eq. (|36|), replacing /(O) and /'(O) by /3(0) and 13(0), which have 
rather complicated explicit forms. We apply these expressions for the two most common 
cases found in Astrophysics: transitions — 1 and 1 — > 2. 

(i) li = 0^lf = l. 

Similarly to what was done in the previous subsection, we consider S{E) for the li = 
— s> // = 1 capture in the case of '^Be(p, 7)"^B" taking a set of binding energies for "^B". 
Results are presented in Fig. |. In this case, there is no initial centrifugal barrier and the 
centrifugal barrier in the final state alone cannot overcome the rising of /|(-E ~ 0) due to 
the pole and the remnant Coulomb. From Eq. (|30|) one can verify that the asymptotics of 
spectral function which goes as (see Eq.p8D, trends to push the peak of the integrand to 
high s- values. This would make hik'^) less sensitive to variation of E, yet it is the third 
derivative that enhances this sensitivity, creating the net effect. 
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By comparing Fig. ^and Fig. ^ one realizes that 5"(0) increases with binding energy for 
both cases. While for = 1 ^ /j = the sign of 5"(0) becomes positive, all = ^ // = 1 
show 5"(0) < 0, suggesting that the centrifugal barrier has a stronger impact in the initial 
state. 

In Fig. ^ the comparison between S{E) and S'(o)(-E') (grey lines) is also done. Whenever 
>lf{e = 0.06, 0.137 MeV) S(o)'(0) > ^'(0); for r]f = lf{e = 0.35 MeV) the two functions 
coincide and for r]f <lf {e = 0.6 MeV) 5(o)'(0) < 5'(0). 

(ii) U = l^lf = 2. 

In Fig. ^ are the S-factors for the /j = 1 — = 2 capture for the reaction ^Be(p, 7)"^B" 
where "^B" is bound by a similar set of binding energies (we chose instead oi e = 0.35 MeV, 
e = 0.09 MeV for which rjf = 2 matching If). As in Fig. |^, 5"(0) > for all cases. Clearly 
here the centrifugal barriers are enough to overcome the influence of the remnant Coulomb 
barrier and the pole ai E = —e. Note that the factor preceding the integral is stronger in 
this case: F{E, 1-^2) = k^^k"^ + r]}K^)II{k^), whereas F{E, 0^1) = k^lKk"^). This factor 
is responsible for the decrease of S(E) as — > 0. 

A subtle feature appears that makes this example contrast with all previously considered 
cases: when the binding energy increases (the pole moves away from the threshold), S'{E) 
decreases. To understand this property, one would have to develop the third derivative 
explicitly. The safe conclusion that can surely be drawn is that the strong centrifugal 
barrier washes out the dependence on binding energy. 

The grey lines in Fig.^ show the 5(0) (E) behavior. Identical conclusions to those obtained 
in (i) can be drawn from the results: when rjf = If S{E) coincides with S(_o){E) (grey circles); 
when If < rjf then 5'(o)'(0) > 5"(0) (solid lines), the opposite happening when If > rjf (dashed 
and dot-dashed lines). 
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C. Charge dependence 



The energy dependence of the S-factors is strongly influenced by the charge of the cap- 
turing nucleus. It is through rji, the remnant Coulomb in the initial state, and rif, the 
Coulomb parameter of the flnal state, that this dependence comes in. In order to il- 
lustrate the relative interplay between these parameters we chose again the typical case 
"W(p,7)"^B", e = 0.137 MeV, where both = // = 1 and = 1 ^ // = 2 tran- 
sitions are considered and the charge of the capturing nucleus is allowed to take the values 
Za = 0,4 and 8 (naturally this last value is unphysical but intends to exaggerate the effect). 

In order to isolate these effects we run three different calculations. First, the S-factors 
resulting from setting Za = 0, such that none of the Coulomb barrier is present. Secondly, 
the results of the calculations for which the remnant initial-state Coulomb barrier is set to 
zero [rji = 0), but keeping rjf. It contains thus only the flnal-state Coulomb barrier. Finally, 
the full calculation for S(E), including both initial and flnal Coulomb barriers. Results are 
plotted in Fig.(P) for a transition 1 — and in Fig.(|^) for a transition 2 — > 1. From these 
flgures we conclude: 

1. the flnal Coulomb barrier increases the S-factor derivative, as a normal penetration 
barrier would; 

2. the larger the charge, the stronger the effect of the flnal-state Coulomb barrier; 

3. the effect of the flnal-state Coulomb barrier is signiflcantly weaker than the initial 
remnant Coulomb barrier; 

4. for all cases, the initial remnant Coulomb barrier signiflcantly decreases the flrst deriva- 
tive, acting as an attractive potential; 

5. the larger the charge, the larger the decrease of 5"(0), such that in most cases it wins 
over all other real barriers making S"(0) negative (note that this happens even for 
li = 2, where all previous cases had 5"(0) > 0). 
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IV. COMPARATIVE BEHAVIOR OF 5^ FACTORS FOR SOME IMPORTANT 

ASTROPHYSICAL PROCESSES 



In this section we compare and explain the difference in the energy behavior of the direct 
capture S-factors for some important astrophysical processes. Since we are only interested 
in the comparative energy behavior, all S-factors are normalized to unity at E = 0. Four 
notorious reactions with either /j = 1 or = 1 are considered: 

a) ^Be(p, 7)^6(2+), e = 0.137 MeV, = ^ // = 1; 

b) i^N(p, 7)1^0(3/2+), e = 0.504 MeV, = 1 ^ = 0; 

c) i'50(p,7)i^F(l/2+), e = 0.105 MeV, 1^ = 1-^1^ = 0; 

d) 20Ne(p,7)2iNa(l/2+), e = 0.0064 MeV, ^^ = 1^/^ = 0. 

Results are presented in Fig.|^. Since /«, If < 1, the centrifugal barriers are not strong enough 
to win over the effect of the singularity a.t E = —e and the initial remnant Coulomb barrier, 
leading to negative slopes for S{0). Within these examples, the largest charge system, also 
corresponds to the lowest binding energy, such that both effects result in a very negative 
derivative for S{0). A significant difference between S"(0) for the *B and ^'^F cases, with 
similar binding energies, proves that the initial Coulomb remnant barrier which is stronger 
for ^^F than for ®B, wins over the initial centrifugal barrier Zj = 1, resulting in a more 
negative derivative S"(0) for ^''F. 

We also consider astrophysical examples for which either /j = 2 or = 2: 

a) ^Be(p,7)8B(2+), e = 0.137 MeV, with 1^ = 2^ If = 1; 

b) i^O(p,7)i^F(5/2+), e = 0.605 MeV, with li = l^lf = 2; 

c) 22Mg(p,7)23Al(5/2+), e = 0.127 MeV, with = 2 ^ = 1. 

The resulting S-factors are presented in Fig.^. For the first two cases the slope is positive. 
For the ®B case, the initial centrifugal barrier and smaller charge result in the highest S"(0). 
The fact that the remnant Coulomb barrier is stronger for ^''F than for ®B and that the initial 
centrifugal barrier is lower, produces a smaller S"(0). For the capture on ^^Na, despite the 
high initial orbital k = 2, the very large initial remnant Coulomb barrier (together with the 
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trace of the nearby singularity) is able to make the slope negative near zero. 

V. SUMMARY 

We investigated the low-energy behavior of the astrophysical S{E) factors for direct 
radiative captures of protons by charged nuclei leading to loosely bound final states. For 
such processes a simple potential model, extending the asymptotic form of the final state to 
r = and assuming a pure Coulomb scattering wave in the initial state are well justified, 
as most of the contribution to the low energy capture to loosely bound states comes from 
large distances. 

We demonstrate that the behavior of the S-factors is governed by six essential ingredients. 
Two act in an attractive sense, creating a negative slope for -S'(0+): the remnant of the initial 
Coulomb barrier (left after extracting the Gamow penetration factor) and the singularity 
at E = —s (where s is the binding energy of the final state). Three act as real penetration 
barriers, producing a positive slope for 5'(0''"): both initial and final centrifugal barriers and 
the final-state Coulomb barrier. The effect of the final centrifugal and Coulomb barriers are 
minor compared with the initial centrifugal and Coulomb. Finally, there is still a photon 
factor /c^^"*"^, which tends to increase the derivative of S{E). We have derived analytical 
expressions for the S-factor in a few typical cases. We have tried to desmystified the idea 
that the energy behavior of the S-factor around threshold is dominated by the pole E = —e. 
We have also shown that by taking only the first term of the Whittaker expansion (final 
state), one cannot reproduce the correct S-factor energy behavior at threshold, except when 
the Coulomb parameters of the final state is integer. Finally, we have not only illustrated 
our findings with a few sets of study cases but applied it to specific examples relevant in 
astrophysics. 

Although in this work, we were mainly concerned with the energy behavior of the as- 
trophysical factors, our approximate equations can reproduce, the absolute values of the 
capture rates within a few percent. The smaller the binding energy, the better accuracy of 
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our equations. The application of Eq. can be extended to energies of a few hundred keV, 
i.e. to the explosive nucleosynthesis energy range, when binding energies are not larger than 
1 MeV. Specifically, for charges = 4 — 10 and binding energies e < 1.0 MeV, S-factors 
obtained from Eq. do not deviate more than a few % from the exact potential model, 
for relative energies up to 1 MeV. However for incident energies > 400 keV and/or higher 
binding energies, a simple potential model becomes less accurate: nuclear interactions in the 
initial state and microscopic effects, such as antisymmetrization, should to be included. 
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APPENDIX A: INTEGRAL REPRESENTATION OF S{E) 



In the Appendix we derive the expression for the S-factor for captures to loosely bound 
states. Our starting point is the definition in Eq. (|^). We simplify this expression, reducing 
it to a first-order integral with an integrand expressed in terms of elementary functions. To 
this end, we use the integral representation for both bound and scattering states, under the 
assumption of the asymptotic approximation for the bound state wave function discussed in 
the main text. 

The Whittaker function has the following integral representation [0]: 

Furthermore, the integral representation of the first term of its asymptotic expansion in 
powers l/ztr (j = in Eq.^), is given by: 

_ — K.r 

riO) 



W"-""' , ,i/o(2Kr) = — 



(JXJ 

y dse-2''"'s'/+^^. (A.2) 





.(0) 



r(// + r// + l) 

The difference between W-r,j:,i^.+i/2{'^ i^r) and the first term W^,^^, i^,^^^2{'^ '^^) integral rep- 
resentations is only in the spectral function. 

The Coulomb regular solution, Fi-{k,r), is given by Eq. (|^) in terms of the confluent 
hypergeometric function. This can be written in integral form as 0: 

iF,{k + l-ir]i, 2h + 2,-2ikr) = . \ ^ . . 

} 1-t 

X dte^"''''[t{l-t)]^^[ -]*''\ (A.3) 



Here, B{li + 1 — irjiji + 1 + irji) is a /3-function expressed in terms of F-functions: 

Substituting Eq. ( |A.1| ) and Eq. ( |A.3D into Eq. (^, allowing for Eq. (^) and changing the 
order of the integrations, we arrive at the following expression for the S-factor: 
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|2 



[r(// + ry/ + i)F [r{2k + 2)] 



X mfk' + vj^') \J\', (A.5) 

We have merged the three integrations into J: 

oo 1 

j-Jds [s{i + s)fnY^/^ J - m^Y"' 





oo 



X J drr^+''+'^+^e~''[''^^+^^^+'''^^~^*^l . (A. 6) 

rO 

We now replace the lower limit of the integration by rg = 0, keeping in mind the application 
to loosely bound states (where the contribution from the interior is not significant). The 
result, after integrating over r, holds: 

S{s) 

oo ^ ^ ^ 1 

j=v Jds [s{i+s)ff / dt[t{i-tt i^r^ 



oo 

with V — J dvv''^~^''f~^^~^^ . Let us define the spectral function of s, as indicated above, 



S{s) and the remaining integral over t as ^7*. The expression for S{s) is: 

= [.(1 + [^^. (A.8) 

The integral over t can be reduced to: 
1 

Jt = l dt[t{l-t)]^^[ 



= (-1) 



t ' [K{l + 2s)+ik{l-2t)fi+h+L+3 
{2 k + 1)1 



X 



{lf + L + li + 2)\ 

H™ }■ 1 — / 1 
lim / dt [t(l - t)^' [^Y^i — ^ ^— — — , (A.9) 



where we introduce m — lf + L + l — li. One can further simplify the integrand by factorizing 
the t-independent terms, such that the remaining integral consists only on the integral 
representation of the hypergeometrical function 2-F1: 
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Ji = (-1) 



(2/. + 1)! 



X lim 



{lf + L + k + 2)\ 
d"" 1 



B{li + l- ir]i,li + l + iT]i) 



Tr^2Fi{k + l-tr]i,2k + 2,2k + 2;y). (A.IO) 



e-o de'" [k(1 + 2s) + e + iA;]2''+2' 

The argument of the hypergeometric Gauss function is y = 1 — Taking into 

account the property: 2F'i{a , P , •y; y) = (1 — 2-^1(7 — Ci,l — P^l'iV) (Eq.9. 131-1 of 

0), and the fact that 2-Fi(«, 0, 7; = 1, we can arrive at a simphfied expression: 

(2/i + l) 



Jt 



-1 



X hm 



{lf + L + k + 2)\ 
d'" 1 



B{li + 1 - ir]i,li + l + ir]i) 



e-o de"^ [{k (1 + 2 s) + e)2 + FJ'.+i 



2?7i arctan , ^ , , 



(A.ll) 



Introducing Eq. ( [A. Ill ) and Eq. ( [A.8| ) into Eq. ( |A.7| ), reduces to 
(2/. + 1)! 



00 



(/j + L + /i + 2) 



B{li + l-ir]iji + l + ir]i) 
1 



X / ds5(s) hm , , , xo ,m, ,1 

J ^ ^ ,^0 de°i [(/«(l + 2s) + e)2 + A;2]'»+i 



^2,7, arctan ^^-,^2s)+e 



(A.12) 



Consequently, replacing Eq. (|A.12|) in Eq. (|A.5|) , S{E) takes the form: 



k2'/+1 '■i 

Vf jJo 



X 



f d"^ 1 

lim dsS(s)- — — — ; : — , , , 

6^07 ^ ^de"! [(/t(l + 2s) + e)2 + P]'«+i 



2 rji arctan , ^ ^ , 

g K(l + 2s) + e 



(A.13) 



where a new constant, putting in evidence energy independent factors, is introduced: 

1 [(2/. + l)!]2 1 



B = A 22''+2'/+4 



TC-. 



(A.14) 



[r(/; + r^f + 1)]2 [{If + L + k + 2)!]2 [r(2 + 2)]2 " 

The constant A was introduced in the beginning of the text by Eq. (|^). 

Similarly, the function S'(o) (E) , resulting from the first term of the Whittaker expansion, 
can be derived. The result holds an identical expression to S{E), the only difference being 
the spectral function over s. Thus, Eq. ( |A.13|) is applicable to S'(o)(-E) when replacing S{s) 
by So{s): 
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(A.15) 



It is important to note that the replacement of the lower integration limit in the radial 
integral of Eq. ( [A .51 ) imposes conditions for the convergence of the integrals over s at cxd: 

1) for S{E), If < 2 + L + li, which is always satisfied; 

2) for S'(o)(-E), rjf < 2 + L + li, which limits the charge of the capturing nucleus. 

Finally, we should mention that, as we focus only on dipole transitions, then |/j — / j| = 1 
and L = 1. Hence, there are two possibilities for m: U — lf = l^m = l and If — U = 
1 — s> m = 3. In each case an evaluation of Eq. ( |A.13| ) can be explicitly performed. This is 



presented and discussed in section |T1. 
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FIG. 1. So{E) for the /j = 0, ^ // = 1 proton capture, for a set of binding energies of the final 
state. Masses and charges are the same as for ''Be(p,7)^B. So{E) is normalized to unity at zero 
energy. 
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FIG. 2. Proton capture h = 1,^ If = for a set of binding energies: S{E) (dark hues) 
versus 5(o)(-E') (grey hues). Masses and charges are the same as for ''Be(p,7)'^B. All S-factors are 
normalized to unity at zero energy. 
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FIG. 3. Proton capture li = 2, ^ If = 1 for a set of binding energies: S{E) (dark lines) versus 
5(0) (-E) (grey lines). For e = 0.09 MeV, S(^q-j{E) is plotted in circles to distinguish from S{E). 
Masses and charges are the same as for ''Be(p,7)^B. All S-factors are normalized to unity at zero 
energy. 
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FIG. 4. Proton capture li = 0, ^ If = 1 for a set of binding energies: S{E) (dark lines) versus 
5(0) (-E') (grey lines). For e = 0.35 MeV, 5(o)(-E') is plotted in circles to distinguish from S{E). 
Masses and charges are the same as for ''Be(p,7)^B. All S-factors are normalized to unity at zero 
energy. 
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FIG. 5. Proton capture /j = 1, — > = 2 for a set of binding energies: S{E) (dark lines) versus 
5(o)(-E') (grey lines). For e = 0.09 MeV, 5(o)(-E') is plotted in circles to distinguish from S{E). 
Masses and charges are the same as for ''Be(p,7)'^B. All S-factors are normalized to unity at zero 
energy. 
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FIG. 6. S-factors for the El-capture ''Be(p, 7)^6: The solid line is the result for S(E) assuming 
Za = 0; the dashed line and dot-dashed correspond to S(E) for = 4 calculated with and 
without the initial Coulomb barrier, correspondingly; the squares and circles correspond to S(E) 
for Za = 8 calculated with and without the initial Coulomb barrier, correspondingly. All S-factors 
are normalized to unity at zero energy. 
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FIG. 7. S-factors for the same reaction as Fig.^, but k = 2, — > If = 1: The sohd hne is the result 
for S(E) assuming Za = 0; the dashed hne and dot-dashed correspond to S(E) for = 4 calculated 
with and without the initial Coulomb barrier, correspondingly; the squares and circles correspond 
to S(E) for Za = 8 calculated with and without the initial Coulomb barrier, correspondingly. All 

S-factors are normalized to unity at zero energy. 
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FIG. 8. S-factors for some astrophysical reactions: '^Be(p, 7)^B(2^), corresponding to a transi- 
tion li = 0^lf = l with binding energy e = 0.137 MeV (solid); ^''N(p, 7)^^0(3/2+), corresponding 
to a transition /j = 1 ^ = with binding energy e = 0.504 MeV (dotted); ^^0(p, 7)^'^F(l/2+), 
corresponding to a transition /j = 1 ^ = with binding energy e = 0.105 MeV (dashed); 
^^Ne(p, 7)^^Na(l/2^), corresponding to a transition = 1 — > = with binding energy e = 0.0064 
MeV (dot-dashed). All the astrophysical factors are normalized to unity at zero energy. 
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FIG. 9. S-factors for some astrophysical reactions: ^Be(p, 7)^B(2+), corresponding to a tran- 
sition li = 2^lf = l and binding energy of the final state e = 0.137 MeV (solid); 
^^0(p, 7)^^F(5/2"'"), corresponding to a transition li = 1 If = 2, and binding energy of the final 
state £ = 0.605 MeV (dashed); ^^Mg(p,7)^3Al(5/2+), corresponding to a transition = 1 ^ // = 2 
and binding energy of the final state e = 0.127 MeV (dot-dashed). All the astrophysical factors 
are normalized to unity at zero energy. 
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